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ABSTRACT 

We present results from monitoring observations of the gravitationally lensed quasar RX J1131 — 1231 
performed with the Chandra X-ray Observatory. The X-ray observations were planned with relatively 
long exposures that allowed a search for energy-dependent microlensing in the soft (0.2-2 keV) and 
hard (2-10 keV) light curves of the images of RX J1131 — 1231. We detect significant microlensing in the 
X-ray light-curves of images A and D, and energy- dependent microlensing of image D. The magnification 
of the soft band appears to be larger than that in the hard band by a factor of ~ 1.3 when image D 
becomes more magnified. This can be explained by the difference between a compact, softer-spectrum 
corona that is producing a more extended, harder spectrum reflection component off the disk. This is 
supported by the evolution of the fluorescent iron line in image D over three consecutive time-averaged 
phases of the light curve. In the first period, a Fe line at E = 6.36^q ies keV is detected (at > 99% 
confidence). In the second period, two Fe lines are detected, one at E — 5.47^0.08 keV (detected at > 
99% confidence) and another at E = 6.02l ; 07 keV (marginally detected at > 90% confidence), and in 
the third period, a broadened Fe line at 6.42+0^5 keV is detected (at > 99% confidence). This evolution 
of the Fe line profile during the microlensing event is consistent with the line distortion expected when 
a caustic passes over the inner disk where the shape of the fluorescent Fe line is distorted by General 
Relativistic and Doppler effects. 



Subject headings: galaxies: active — quasars: individual (RX J 1131-1231) 
gravitational lensing 



X-rays: galaxies 



1. INTRODUCTION 

Some of the most fascinating properties of quasars arise 
from the strong gravitational field that dominates over all 
other forces near the supermassive black hole. Under the 
effects of strong gravity and rapid rotation, the energies 
and kinematics of X-ray photons emitted at a few gravita- 
tional radii are modified significantly. Observations of the 
X-ray spectra of quasars thus have the potential of test- 
ing the theory of General Relativity in a region of strong 
gravity, constraining the structure of the emitting region 
and providing estimates of the spin of the black hole (e.g., 
Fabian et al 1989; Laor 1991). 

A promising technique for indirectly mapping the ac- 
cretion disk around the black hole, and thus constraining 
the properties of the black hole, relies on monitoring varia- 
tions in the profiles of Fe fluorescence lines originating from 
the inner parts of the disk (e.g., Fabian et al 1989; Laor 
1991). The iron fluorescence line is produced by X-rays 
from a hot corona irradiating the accretion disk. The line 
typically contains a narrow and a broad component. The 
narrow component is thought to originate from material in 
the broad-line region and/or a molecular torus. The broad 
component most likely originates from accreting material 
closer to the black hole. The iron line method has been 
most successful in studies of Galactic black-hole binaries 
where the X-rays are mostly thermal emission from the 
accretion disk (e.g., Done, Gierlihski, & Kubota 2007 and 
references therein). For the cooler, quasar disks, the X- 
rays are dominated by the non-thermal emission from the 



corona rather than direct emission from the disk, making 
it more difficult to understand and model. 

The broad line profiles are determined by a combina- 
tion of the Doppler, special and general relativistic (GR) 
effects (e.g., Fabian et al 1989; Laor 1991), while their vari- 
ations are caused by fluctuations in the X-ray continuum 
that drives the line emission (Young & Reynolds 2000). 
Modeling of the distorted iron line can provide informa- 
tion about the disk geometry, inclination and the spin of 
the black hole (e.g., Fabian et al. 2000; Reynolds & Nowak 
2003). The accuracy of the Fe line method for constraining 
the spin of a black hole from the broadened red wing of 
the Fe line profile has been questioned recently (e.g., Miller 
et al 2008; Sim et al. 2010). In particular, simulations by 
Sim et al. (2010) show that the Fe line can develop a weak, 
red-skewed line wing as a result of Compton scattering in 
an outflow. Separating the broadening of the Fe line by 
GR effects near the inner edge of the disk from that pro- 
duced by scattering in a wind may be difficult, especially 
with the limitations of low resolution X-ray spectra. 

An independent method of probing the structure of 
quasars relies on gravitational microlensing (e.g., Paczyn- 
ski 1986; Wambsganss 1990). Stars in the foreground lens 
galaxy produce complex microlensing magnification maps 
that pass over the quasar. The characteristic size of a mag- 
nification caustic produced by a single star is referred to 
as the Einstein radius. The strength of the magnification 
depends on the relative size of the emitting region and the 
Einstein radius, with more compact sources having higher 
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microlensing variability amplitudes. 

Microlensing observations constrain the sizes of accre- 
tion disks at optical wavelengths (Pooley et al. 2007; An- 
guita et al. 2008; Morgan et al. 2010; Jimenez- Vicente 
et a. 2012), the sizes of quasar coronae at X-ray ener- 
gies (Chartas et al. 2002, 2007, 2009; Pooley et al. 2007; 
Morgan et al., 2008; Dai et al. 2010; Blackburne et al. 
2011b) and the dependence of the accretion disk size with 
wavelength (Poindexter et al. 2008; Bate et al. 2008; 
Eigenbrod et al. 2008; Floyd et al. 2009; Morgan et al. 
2010; Blackburne et al. 2011a; Mosquera et al. 2009, 
2011b). A study of microlensing variability in a sample 
of eleven microlensed quasars performed by Morgan et al. 
(2010) showed that the radius of the accretion disk scales 
as and the sizes estimated from microlensing are 

larger than predicted by thin disk theory for the typical 
Eddington ratios, (ieoi/^Edd ~ 1/3), of quasars. 

We recently presented the results of a six epoch moni- 
toring campaign of the z s =0.658 quasar RX J1131 — 1231 
(Chartas et al. 2009; Dai et al. 2010). In these studies we 
found significant flux variability in all four images and that 
the saddle point image A was strongly microlensed (also 
see Blackburne et al. 2006; Pooley et al. 2007). In Dai et 
al. (2010) we modeled the X-ray and rest-frame UV light 
curves of RX J1131 1231 in detail using the technique in- 
troduced in Kochanck (2004), and found that the X-ray 
emission in RX J1131 — 1231 comes from a region smaller 
than ~ 10r 9 while the observed optical emission region has 
a size of ~ 70r g , where r g = GMbh/c 2 is the gravitational 
radius of the black hole and RX J1131 — 1231 is estimated 
to have Mbh = 1-3 x 1O 8 M based on the Hp line width 
(Peng et a. 2006). 

Here we present the results from the continuation 
of our X-ray and optical monitoring campaigns for 
RX Jl 131 — 1231. The X-ray exposure times in the new 
epochs have been increased to provide sufficient signal-to- 
noise spectra to search for X-ray energy dependent mi- 
crolensing. In §2 we present the analysis of the X-ray ob- 
servations, in §3.1 we discuss the detection of microlensing 
events in the light curves of images A and D, in §3.2 we 
present spectra possibly showing the evolution of the Fe 
line as a caustic crosses the black hole, and in §3.3 we 
present constraints on quasar structure assuming the de- 
tected Fe line energies are correct. Throughout this paper 
we adopt a flat A cosmology with H = 70 km s _1 Mpc -1 , 
n A = 0.7, and M = 0.3. 

2. X-RAY OBSERVATIONS AND DATA ANALYSIS 

RX J1131-1231 was observed with the Advanced CCD 
Imaging Spectrometer (ACIS; Garmire 2003) on board the 
Chandra X-ray Observatory (hereafter Chandra) 29 times 
between 2004 April 12 and 2011 February 25. Results from 
the analysis of the first six observations of RX J1131 — 1231 
were presented in Chartas et al. (2009) and Dai et al. 
(2010). Here we describe the data analysis for the remain- 
ing 23 observations, but we will use the results from all 29 
observations in our microlensing analysis. 

The Chandra observations of RX J1131 — 1231 were an- 
alyzed using the standard software CIAO 4.3 provided by 
the Chandra X-ray Center (CXC). A log of the observa- 
tions that includes observation dates, observational iden- 
tification numbers, exposure times, ACIS frame time and 



the observed 0.2-10 keV counts is presented in Table 1. 
We used standard CXC threads to screen the data for sta- 
tus, grade, and time intervals of acceptable aspect solution 
and background levels. To improve the spatial resolution 
we employed the sub-pixel resolution technique developed 
by Tsunemi et al. (2001). 

We estimate the X-ray counts of the images in the total 
(0.2-10 keV), soft (0.2-2 keV) and hard (2-10 keV) bands. 
The total counts were corrected for pile-up for images A 
and B, and pile-up is unimportant for images C and D. The 
hard and soft bands for images A and B were not corrected 
for pile-up because of the difficulty in estimating energy- 
dependent pile-up corrections at the observed count rates. 
Pile-up is an instrumental effect that occurs when two or 
more X-ray photons strike individual or neighboring CCD 
pixels within one frame time and can cause spectral and 
point-spread-function (PSF) distortions, grade migration, 
and an apparent hardening of the energy spectra. We 
used the forward spectral-fitting tool LYNX (Chartas et 
al. 2000) to estimate the fraction of events lost due to the 
pile-up effect. These corrections affect image A (up to 53% 
in epoch 23) and B (up to 22% in epoch 18). Pooley et al. 
(2012) in their analysis of the epochs from 2004 April 12 
to 2008 April 11 did not include these corrections, which 
may explain some of the epoch-dependent results in their 
analysis. 

We fit the lensed images simultaneously using the rela- 
tive astrometry derived from HST images (Morgan et al. 
2006) and MARX models (Wise et al. 1997) of the PSF. 
We fit the data using 0';0246 bins (compared to the 0'f491 
ACIS pixel scale) by minimizing the C-statistic (Cash 
1979) between the observed and model images. The PSF 
fitting method provided the best fit value and errors for 
the number of counts in each image, Ai jPS f and cr i;Psf , re- 
spectively. We define as A^ABCD, ps f = J2i ^i.psf to be the 
total counts of all images from the PSF fitting method. 

While the PSF models should correctly recover the rel- 
ative fluxes, errors in the PSF models may bias the total 
fluxes. We corrected for this by renormalizing the total 
counts to match the total counts found in a 5" radius re- 
gion centered on the lens with a background correction 
based on an annulus from 7'.'5 to 50" around the lens. 
Specifically, the counts from each image as determined 
from the PSF fitting method were renormalizcd to N{ { = 

RxiV i;PS f, where R = A A BCD,a P /A A BCD,PSF and iVABCD, ap 
is the total counts from the aperture method. Table 2 lists 
the exposure times and the N[ f counts of images C and 
D in the soft and hard bands. 

Figure 1 shows the stacked image of all observations 
of RX J1131-1231 listed in Table 1, with a total effec- 
tive exposure time of ~ 333.6 ks. To reduce background 
contamination and to enhance possible soft extended X- 
ray emission we filtered the data to include only photons 
with energies from 0.4 and 3.0 keV. The image was binned 
with a bin size of 0.5 arcsec and adaptively smoothed with 
the CSMOOTH tool developed by Ebeling et al. (2000). 
CSMOOTH smooths a two-dimensional image with a cir- 
cular Gaussian kernal of varying radius. We find no emis- 
sion from a possible cluster of galaxies in the field with 
a 3 g upper limit on the 0.2-2 keV luminosity of 6.1 x 
10 41 ergs s _1 , assuming a cluster with a temperature of 
1.5 keV at the lens redshift, an abundance of 0.3 solar and 
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an extraction radius of 30 arcsec. 

The stacked image of RX J1131 — 1231 does not show 
any additional lenscd images. To obtain quantitative lim- 
its on possible additional images located at the center 
of RX J1131-1231 we extracted the 0.2-10 keV counts 
within a 0.4 arcsec circle centered on the four images and 
the lens. The backgrounds in these extraction regions are 
dominated by the contamination from the images rather 
than the instrumental or cosmic backgrounds. To esti- 
mate the fractional contamination per image we used the 
0.2-10 keV counts in 0.4 arcsec circular apertures placed 
1 and 2 arcsec North of image B in RX J1131-1231. Be- 
cause the summed image averages over many spacecraft 
orientations, the averaged PSFs are largely free of angular 
structure. The detected 0.2-10 keV counts in the 0.4 arc- 
sec aperture centered on the lens are consistent with the 
estimated background and contamination from the bright 
images. Combining all the data we obtain a 3<r upper 
limit on the 0.2-10 keV flux of any central, odd image of 
4.7 x 10~ 15 ergs s _1 . This limit corresponds to a 3<r limit 
on the flux ratio of a factor of ~ 135 relative to the aver- 
age flux of image C during period 2. While this is a strong 
limit on the flux of a central image, it is not yet strong 
enough to place interesting constraints on the central den- 
sity of the lens (see Keeton 2005). 

In Figure 2 we show the 0.2-10 keV light-curves of im- 
ages A, B, C and D. Images A and D show significant 
uncorrelated flux variability suggesting that these images 
are also significantly affected by microlensing. The light- 
curves have been shifted to account for the estimated time 
delays between the images of RX J1131 — 1231 (Chantry 
et al. 2010). The lower amplitude variability seen in im- 
ages B and C is strongly correlated, suggesting that these 
images are less affected by microlensing. 

In Figure 3 we show the observed X-ray and optical flux 
ratios of D/B, A/B, and C/B. The i?-band monitoring 
data shown in Figure 3 are obtained with the SMARTS 
1.3m and ANDICAM (Depoy et a. 2003). They are re- 
duced using the standard procedures and modeled as de- 
scribed in Kochanek et al. (2006). The time delays be- 
tween images are considerably shorter than the time scales 
of the two major magnification events seen in images A 
and D and the intrinsic brightness of the source is fairly 
constant over the observational period, as seen in the light- 
curves of images C and B. We therefore do not expect the 
observed flux ratios to be significantly affected by intrin- 
sic variability over the period of the strong magnification 
events in images A and D. 

The total X-ray band flux ratio C/B is relatively con- 
stant at 0.29 ± 0.03, supporting our hypothesis that the 
microlensing magnifications of images B and C are rela- 
tively constant during our new observations. It also indi- 
cates that the intrinsic flux of the quasar does not vary by 
more than 28% during the observing period. Specifically, 
the ratios of the standard deviations to the mean total 
band X-ray count rates of images A, B, C and D are 0.58, 
0.28, 0.24, and 0.75, respectively. The large flux variations 
detected in images A and D are uncorrelated with those 
of C and B and are therefore not intrinsic to the quasar 
but must be the result of microlensing. 

Image A is significantly and image B is mildly piled- 
up, and, as a result, their spectra are expected to appear 



flatter in energy and distorted. We therefore focus the 
spectral analysis on images C and D. We fit the Chandra 
spectra of images C and D of RX J1131 — 1231 for each 
epoch with a model that consists of a power law with neu- 
tral intrinsic absorption at z s — 0.658. We also included 
Galactic absorption due to neutral gas with a column den- 
sity of N H = 3.6 x 10 20 cm~ 2 (Dickey & Lockman 1990). 
The results of these fits are presented in Table 3. As shown 
in Table 3, the values of the photon indices and X-ray 
fluxes of image C do not show any significant variations. 
Significant spectral variability is detected in image D that 
is likely associated with energy-dependent microlensing. 
The energy dependent microlensing interpretation is dis- 
cussed in more detail in the section §3. The spectral fits 
also indicate that there is no significant absorption at the 
lens or source redshifts. 

3. DISCUSSION 

3.1. Microlensing events RX J1131-1231 

Significant deviations of the flux ratios from a constant 
value over time are suggestive of microlensing since the in- 
trinsic variability of the quasar cancels in a flux ratio up to 
our ability to correct for time delays. We observe signifi- 
cant variability in the ratios D/B, D/C, A/B and A/C but 
not in C/B suggesting that only images A and D show sig- 
nificant microlensing activity during our new observations 
(e.g., see Figure 3). Specifically, we find that the total 
band flux ratio A/B deviates significantly from a constant 
between days 4800 and 5600, a period of about 800 days. 
The total flux ratio D/B deviates from a constant for a 
period of about 1,100 days. 

Microlensing variability occurs when the complex pat- 
tern of caustics produced by stars in the lens moves across 
the source plane. The characteristic scale of these patterns 
is the Einstein radius of 
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where (M) is the mean mass of the lensing stars, D is the 
angular diameter distance, and the subscripts L, S, and O 
refer to the lens, source, and observer, respectively. Emis- 
sion regions with sizes significantly larger than Re will be 
affected little by microlensing, whereas emission regions 
significantly smaller than i?E can be strongly magnified. 
Due to the combined motions of the observer, lens, source, 
and stars in the lens galaxy, the source moves relative to 
the magnification patterns at an effective velocity of order 
700 km s' 1 for RX J1131-1231 (Mosquera & Kochanek 
2011a). The time scale to sample a significantly different 
microlensing region is the Einstein crossing time of t& = 
Re/v c ~ 11 years, which is roughly the temporal baseline 
of the X-ray data. Microlensing variability can be seen on 
much shorter time scales corresponding to the time scale 
needed to move by the size of the source, t§ = Rg/v c . For 
the optical (R = 1.3 x 10 15 cm) and X-ray (R x = 2.3 
x 10 14 cm) sizes we found in Dai et al. (2010), the source 
crossing times are of the order to = Ro/ v c — 0.6 years 
and ix = Rx/v c — 0.1 years. 

Figure 3 shows these time scales in comparison to the 
changes in the flux ratios, and it is striking that the ap- 
parent time scales in the X-ray and optical A/B and D/B 
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flux ratios closely match the predictions from the earlier 
size estimates. For comparison, the estimated black hole 
mass of 6 x 10 7 M (Peng et al. 2006) corresponds to a 
Schwarzschild radius of Rs = 1.8 x 10 13 cm. Microlensing 
magnification patterns have, however, strong spatial cor- 
relations with dense networks of high magnification caus- 
tic regions separated by broad valleys with little activity. 
Thus, an image showing strong microlensing variability is 
likely to continue to do so for a time scale of order t^, 
while one without such variability is also likely to con- 
tinue in such a state. Previous studies of RX J1131 — 1231 
suggest that images A and D are "active" and images C 
and B are not. 

3.2. Energy Dependent Microlensing in RX J1131-1231 

Figure 4 shows the soft (0.2-2 kcV) and hard (2-10 kcV) 
flux ratios D/B and D/C. The total, soft, and hard band 
flux ratios D/C and D/B are almost constant with time 
between days ~ 4000 — 4500. After this period the flux ra- 
tios increase, indicating the onset of significant microlens- 
ing magnification of image D. The soft and hard band flux 
ratios begin to differ at a modified Julian date of ~ 4,700 
days, just after an increase in the microlensing magnifica- 
tion of image D. Specifically, the soft band flux ratios D/B 
and D/C are systematically larger than the hard band ra- 
tios by a constant factor of 1.3 ± 0.1 for a period of about 
600 days. We therefore associate these differences of the 
soft and hard band flux ratios D/B and D/C with an en- 
ergy dependent microlensing event in image D. Image A 
is significantly affected by pile-up and any estimate of the 
soft and hard band fluxes of this image is very uncertain. 
We have therefore not attempted to correct the soft and 
hard band fluxes of image A, but only its total band flux. 

It is possible to explain the energy dependent behavior 
of image D as a caustic approaching and traversing the 
accretion disk. The X-ray emission is sufficiently compact 
that its microlensing variability is probably dominated by 
the effects of individual caustics (although there will be 
exceptions). We assume the X-ray spectrum is a combina- 
tion of a direct X-ray component from a hot corona that is 
centered on the black hole and a more extended reflected 
component from the accretion disk. As shown in Figure 5, 
adapted from Reynolds et al. (1998), the reflected X-ray 
component is harder than the coronal component. The 
ratio of rcflected-to-direct emission is larger in the hard 
band than in the the soft band, 

(7^0 >{4^) > (2) 

\ J Direct J Hard V J Direct J Soft 

where /roA is the flux of reflected emission off the accre- 
tion disk and /Direct is the flux of direct emission from the 
hot corona. 

When the central source is far away from a caustic, 
the two components have similar magnifications and the 
observed spectrum is dominated by direct emission from 
the corona with a smaller contribution from reflected disk 
emission. As the caustic approaches the center of the black 
hole, the more extended reflected X-ray component is mag- 
nified first. The X-ray spectrum at this point is still domi- 
nated by the direct component and the ratio of soft-to-hard 
band is not expected to change significantly. As the source 
approaches the caustic, the more compact corona becomes 



strongly magnified. Since the coronal spectrum is softer 
than the reflected, the soft band will be more magnified 
than the hard band as the source crosses and passes the 
caustic. 

We can use the evolution of the X-ray spectrum of image 
D to test our explanation. To improve the S/N we stacked 
the X-ray spectra of images C and D in three periods. 
Period 1 includes epochs 7 though 16, period 2 includes 
epochs 17 through 22, and period 3 includes epochs 23 
through 29. In Figure 6 we overlay the stacked spectra of 
images C and D for the three periods to show the spec- 
tral evolution of the continuum and the lines. The stacked 
spectra of image D are also displayed in separate panels 
in Figure 7 to clearly show the iron lines. In Figure 8 we 
show the 68%, 90% and 99% confidence contours of the 
flux normalizations of the Fe lines detected in images C 
and D versus their energies during the three periods. 

During period lthe D/C flux ratios remain constant 
given the uncertainties, indicating that the microlensing 
magnifications are roughly constant. A significant Fe K 
line is detected at an energy of 6.36lg!i6 keV that is unre- 
solved in width. The 0.2-10 keV flux in the entire spectrum 
of image D for period 1 is 2.1 x 10~ 13 erg s" 1 cm" 2 , the 
flux in the Fe line is 4.2 x 10~ 15 erg s _1 cm~ 2 and its 
rest-frame equivalent width is We = 398^8 eV. During 
period 2, the D/C flux ratios begin to increase, indicating 
a change in the magnification of image D, but no energy 
dependence is detected in the flux ratios. A significant Fe 
K line is detected at an energy of 5.47+ ' 08 keV (unresolved 
in width) and a second line is marginally detected at an 
energy of 6.02tg;°y keV (unresolved in width). The 0.2- 
10 keV flux in the entire spectrum of image D for period 
2 is 2.7 x 10~ 13 erg s _1 cm~ 2 and the fluxes (rest-frame 
equivalent widths) in the 5.48 keV and 6.02 keV lines are 
3.4 x 10~ 15 erg s" 1 cm" 2 (W E = 210±^ eV) and 3.0 x 
10~ 15 erg s" 1 cm- 2 (WE = 200±i|^ eV), respectively. Dur- 
ing period 3, the D/C flux ratios are significantly enhanced 
and the soft band flux ratios are larger than the hard band 
ratios by a factor of 1.3 ± 0.1. A significant broad Fe K 
line is detected at an energy of 6.42+°; 15 keV with a width 
of 220^ eV. The 0.2-10 keV flux in the entire spectrum 
of image D for period 3 is 7.5 x 10~ 13 erg s _1 cm~ 2 and 
the flux (rest-frame equivalent width) in the Fe line is 1.2 x 
10~ 14 erg s- 1 cm- 2 (VF E = 290±^g eV). The spectrum of 
image C is essentially unchanged over the three periods, as 
shown in Figure 6. In Figure 9 we show the stacked spec- 
trum of image C from combining all epochs. A significant 
Fe K line is detected at an energy of 6.371°;°' keV. The 
average 0.2-10 keV flux in the combined spectrum of im- 
age C over the three periods is 6.1 x 10~ 13 erg s _1 cm~ 2 
and the flux (rest-frame equivalent width) in the Fe line is 
4.5x 10" 15 erg s" 1 c m - 2 (VF E = U0±f eV). 

The X-ray data imply that the ratio of the total flux 
(corona dominated) to Fe line flux increases from Period 
2 to 3, which would imply that the spectrum of D would 
become softer from Period 2 to 3 since the coronal spec- 
trum is softer than the reflected one. In Figure 10 we show 
the ratios of the spectra and the best-fit absorbed power- 
law plus line models of images D and C. Specifically, we 
plot the observed and modeled ratios D(E)i/C(E), where, 
Di(E) is the flux density for image D in period i and C(E) 
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is the flux density for image C in all three periods. The 
three periods are shown in Figure 4. We detect a signif- 
icant softening of the spectrum during the onset of the 
strong microlensing event in image D that occurs between 
Periods 2 and 3. We note that the softening of the spec- 
trum of an image during a microlensing event is opposite 
to what we observed in Q2237+0305 (Chen et al. 2011). 

3.3. General and Special Relativistic Effects 

Assuming both detections of Fe line components in 
image D during period 2 are correct, we can interpret 
them as the result of distortions of the fluorescent Fe line 
due to general relativistic (GR) and special relativistic 
(SR) Dopplcr effects. Detailed models are likely unwar- 
ranted given the present data, so we will focus on a semi- 
quantitative outline of what we have observed. Assuming 
that the Fe Ka emission arises from the disk, the energy 
shift of the midpoint of the two lines can be used to con- 
strain the radius of the emitting material or equivalently 
the distance of the caustic from the black hole. The split- 
ting of the lines is likely the result of the Doppler effect 
arising from an approaching and receding part of the disk. 
Both lines are expected to be redshifted by the GR ef- 
fect but also for large velocities special relativity predicts 
a redshift of both lines. 

The gravitational redshift of the Fe line is 



1 - 2 



^cm 



(3) 



where r g is the gravitational radius, r cm is the radius of 
the line emitting region and Eq — 6.4 keV. The redshift 
due to special relativistic effects (relativistic Dopplcr) of 
the midpoint of the two lines is 



AEsr — Eq 
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1 



(4) 



Here AEsr only includes the redshift of the midpoint of 
the two lines and not their splitting. We have assumed 
Keplerian velocities with (wfe/c) 2 = r g /r em for simplicity. 
The predicted total shift of the midpoint between the two 
lines from E = 6.4 keV is AE To tai = AE GR + AEsr. 
We observe a total energy shift AExotai = 0.65 keV. 
This energy shift implies (based on eqs. [3] and [4]) that 
r em <~ 15r ff and the corresponding Keplerian velocity at 
this radius is 0.26c. The estimated value of the Keplerian 
velocity, v^, combined with the observed separation of the 
two lines provides a constraint on the inclination angle, 
Unci , of the disk. Specifically, the energy separation of the 
two lines is 



cos 9 



(5) 



where, 7 is the Lorentz factor and ii nc i = 90° — 6. The 
observed line separation of 0.55 ± 0.10 keV implies (from 
cq. [5]) an inclination angle ii nc i = 10^2 degrees. 

The direction of motion of the caustic with respect to 
the projected disk axis will affect the relative amplitude of 
the lines. This may explain why the Fe line from the re- 
ceding material at E = 5.5 keV is brighter than the line at 



E = 6 keV from the approaching material. The redshifted 
lines are narrow and not broad as typically observed in 
Scyfcrt galaxies. This is likely the result of the caustic 
selectively magnifying emission within a narrow radius as 
it approaches the black hole. 

As we showed in Figure 6 the continuum and the line 
emission are relatively unchanged between epochs in the 
non-microlensed image C. These spectra confirm that the 
quasar is not significantly showing any intrinsic spectral 
variability between observing periods. The detected shift 
of the Fe line in image D of RX J1131 — 1231 is therefore 
not the result of an intrinsic phenomenon such as a rotat- 
ing hotspot. This reinforces our conclusion that the shift 
in the energy of the iron line in the microlensed image D 
during period 2 is the result of selective magnification of 
material closer in to the black hole as the caustic crosses 
the accretion disk. 

We have observed time-dependent X-ray microlensing 
in all seven systems examined in our Chandra monitoring 
program (Chen et al. 2012). We see energy dependent 
effects in five systems, although in all cases the differences 
are modest in the sense that the differences between the 
hard and soft bands are far smaller than the differences 
between the X-ray and optical microlensing. Thus, while 
the hard and soft X-ray emission regions may have differ- 
ent average sizes, those differences are not large. Similarly, 
all our analyses to date indicate that the X-ray emission 
region has a size of order 10r g . Our upper limits on the X- 
ray sizes are robust, however, our lower limits are currently 
not well constrained. An improvement on the lower lim- 
its requires stronger limits on the maximum possible mi- 
crolensing variability amplitude, information that is easily 
obtained with additional monitoring of these objects. The 
observed changes in the Fe Ka line structure implies that 
the Fe emission region of RX J1131 — 1231 is compact, re- 
gardless of the details of our simple quantitative model. 
For the Fe line study, the challenge is simply to obtain 
enough X-ray photons to fit the Fe line with more realistic 
models that can provide constraints on the structure of the 
disk and the properties of the black hole. A strategy than 
can address both the lower bound and the Fe line prob- 
lems is to observe lenses similar to RX J1131 1231 that 
have a high probability of showing Fe line microlensing 
using a denser time sampling of epochs that are individu- 
ally long enough to measure the hard and soft band fluxes 
accurately. The denser time sampling would solve the min- 
imum source size problem. Co-adding the epochs would 
then provide a study of the Fe Ka line at the cadences we 
presently achieve for the continuum flux ratios. 
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TABLE 1 

Log of Observations of Quasar RX J1131 — 1231 









Chandra 


Exposure 












Epoch 


Observation 


JD a 


Observation 


Time 




N A C 


JV B C 




N D c 




Date 


(days) 


ID 


(ks) 


00 


counts 


counts 


counts 


counts 


1 


2004 April 12 


3108 


4814 


10.0 


3.14 


425+^ 2 


295OI5J 


839 




2 


2006 March 10 


3805 


6913 


4.9 


0.741 


nnn+20 

393±20 


624j; 2 5 


2041H 


1031 IS 


3 


2006 March 15 


3810 


6912 


4.4 


0.741 


ooi +20 

381 + 20 


6I6I25 


2331H 


93l 10 


4 


2006 April 12 


3838 


6914 


4.9 


0.741 


A 1 Q+20 

413 -20 


507 -23 


1461 - 


1 01 +12 

131-12 


5 


2006 November 10 


4050 


6915 


4.8 


0.741 


3708lgJ 


villas 


36711^ 


155111 


6 


2006 November 13 


4053 


6916 


4.8 


0.741 


3833lg2 


1618l 40 


A n r- + 20 
415 -20 


H5IH 


7 


2006 December 17 


4087 


7786 


4.88 


0.841 


or A i +102 

3541I 9 g 


1443T 5S 


4171 28 


1 1 -,+ 15 


8 


2007 January 01 


4102 


7785 


4.70 


0.441 


2305^74 


1082l 4 g 


3121 25 


108ll| 


9 


2007 February 13 


4145 


7787 


4.71 


0.441 


, c1 +79 

2451 + 78 


IH6I49 


sollal 


1691 I7 


10 


2007 February 18 


4150 


7788 


4.43 


0.441 


2232^73 


9501S 


252l 22 


., ., „ + 15 
11511° 


1 1 


2007 April 16 


4207 


7789 


4.71 


0.441 


o ■"> o o + 74 

23281 74 


1 nnn + 54 

1202l 53 


3441 2 6 


qq+14 

yy -12 


12 


2007 April 25 


4216 


7790 


4.70 


0.441 


2043_ 70 


i r\c o + 50 

10631 48 


3631 ae 


142 -16 


13 


2007 June 04 


4256 


7791 


4.66 


0.441 




14801™ 


373I27 




14 


2007 June 11 


4263 


7792 


4.68 


0.441 


22541^2 


„ , „ „ 4- fi 8 

1466lJ° 


+26 

337l 2 ° 


, 4.1 a 
1291 Is 


15 


2007 July 24 


4306 


7793 


4.67 


0.441 


1958l 68 


1 OO A +57 


353I25 


811n 


16 


2007 July 30 


4312 


7794 


4.67 


0.441 


2725±to 


-■ r-, A A 4-fiQ 

18441^ 


4961^ 


- ^.-,4-1 ^ 

100IH 


17 


2008 March 16 


4542 


9180 


14.32 


0.741 


5557I 117 


43471 104 


13371^ 


351124 


18 


2008 April 13 


4570 


9181 


14.35 


0.741 


n i on + 147 

8199I 147 


r-^ r- A +118 


14531^ 


377l 24 


19 


2008 April 23 


4580 


9237 


14.31 


0.741 


6786+ 130 


4927^^ 


1 i-lTi-l + 5 1 

1279_ 51 


2321^ 


20 


2008 June 01 


4619 


9238 


14.24 


0.741 


4647 j;^® 


3252l 88 


0-70+43 

878+ 41 


463l 2 ? 


21 


2008 July 05 


4653 


9239 


14.28 


0.741 


5587+;i! 

auo ' —118 


3584+i^ 


1001+" 


635 + 32 


22 


2008 November 11 


4782 


9240 


14.30 


0.741 


113 


3085 + |g 


885ll 2 


4881™ 


23 


2009 November 28 


5164 


11540 


27.52 


0.741 


360241^2 


73571-2? 


24201^1 


3827i§ 1 ' 


24 


2010 February 09 


5237 


11541 


25.62 


0.741 


26850^290 


5814+H, 


205915*; 


24371™ 


25 


2010 April 17 


5304 


11542 


25.67 


0.741 


20935±2^ 


61241-° 


1962+g 2 


28131™ 


26 


2010 June 25 


5373 


11543 


24.62 


0.741 


1852lt 2 22 ° 


54451^ 


1522l|* 


1487l|l 


27 


2010 November 11 


5512 


11544 


25.56 


0.741 


270771^ 


6316lll| 


182ll 64 


1730+|= 


28 


2011 January 21 


5583 


11545 


24.62 


0.741 


5689ti^ 


3175_ 123 


10081^ 


982lt 6 


29 


2011 February 2 


5618 


12833 


13.61 


0.741 


64121 1 " 


4375ll 28 


H6III2 


8211" 



a Julian Date — 245000. b ACIS frame-time. c Background-subtracted source counts for events with energies in the 0.2-10 keV band. The counts for 
images A and B of RX J1131 — 1231 are corrected for pile-up. 
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TABLE 2 

Soft (0.2-2 keV) and Hard (2-10 kcV) counts of images C and D 



Observation 
Date 


Soft 
N c 
counts 


Soft 
counts 


Hard 

N c 
counts 


Hard 

N-d 
counts 


2004 April 12 


6621* 


1721H 


— 1 --4-I3 

1 1 1 -13 


3911 


2006 March 10 




__j.q 


-~ + 7 
5OI7 


2411 


2006 March 15 




7211 


6i±t 


2211 


2006 April 12 


10611° 


1031H 


4011 


31±| 


2006 November 10 


28et]_l 


117111 


81±| 


3811 


2006 November 13 


320+H 


981 10 


961 10 


1 T + 4 

17l 4 


2006 December 17 


3381 If 


lOOlg 


90111 


20l| 


2007 January 01 


260ll^ 




55l 8 8 


20l| 


2007 February 13 


_ . _ _i_i7 

249l 16 


144 -12 


58lg 


2911 


2007 February 18 


197ll| 


9411 


5811 


22l° 


2007 April 16 


270II7 




79 1 g° 


16l* 


2007 April 25 


272!!^ 


1 1 ,-+11 


90-10 


27 +6 


2007 June 04 


300_Jg 


,,,4.9 
83ls 


75lg 


2611 


2007 June 11 


28211? 


t m +10 

IOII10 


™ + 9 

5911 


2911 


2007 July 24 


28511? 


^ a + 8 


-,,4-q 

7211 


-1 ^ + 5 

1811 


2007 July 30 


392I2I 


f, A -l-Q 

84lg 


106llj 


_ _-|_4 


2008 March 16 


IO58I34 


2561H 


2851H 


97Tg 


2008 April 13 


, , ~„ + 36 

1187l|° 


30811 J 


29211? 


7611 


2008 April 23 


972lg 


1691H 


, , 1 _+20 

3 171 20 


^1- -1-8 

6511 


2008 June 01 


6721^ 


364^ 


217+11 


1051H 


2008 July 05 


811+ 29 

-30 


492 +21 


202 + l^ 


151 + 12 
±iJ i_ 12 


2008 November 11 


652lf 6 


000 + 20 
JSS_ 2Q 


24511? 


1071H 


2009 November 28 


1952lf 6 


322ll" 


6291* 


866111 


2010 February 09 


1588l« 4 9 


2017111 


5581H 


5231H 


2010 April 17 


15411" 


22241*1 


4841* 


681IH 


2010 June 25 


12131H 


11581H 


3681™ 


00 +20 

d»»_20 


2010 November 11 


14601J1 


1359l|l 


453111 


459111 


2011 January 21 


767111 


7251=1 


251 + 1? 


268l 2 ! 


2011 February 2 


922111 


619125 


2501- 


209111 
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TABLE 3 

Spectral Analysis of Image C of RX J1131-1231 



Observation 




r 


7V H (^=0.658) a 


x7^ 




Soft Flux 


Hard Flux 


Date 






(10 22 cm -2 ) 






(10~ 13 ergs 


cm -2 s~ 1 ) 


2004 April 12 


1 


70 

■ lu -0.11 


_ ^+n oe 
O-OIo.o 


54/56 


0.57 


3.14 


3.79 


2006 March 10 


1 


■70+0.29 
■ ' — 0.26 


,-. „-|_n 1 7 

u - u -o.o 


14/15 


0.52 


1.26 


2.31 


2006 March 15 


1 


OU -0.22 


n „-un so 

u - u -o.o 


18.1/16 


0.36 


1.84 


2.95 


2006 April 12 


1 


_ 1 _|_o 51 
■ bi ~-0.35 


^ ^ + 35 

u - u -o.o 


7.4/9 


0.60 


1.22 


1.53 


2006 November 10 


1 


■ y,:i -0.19 


^ ^ + 09 

0-0Io.o 


34/27 


0.16 


3.36 


2.84 


2006 November 13 


1 


qq+0.24 
' yy -0.21 


n i +0.2 

O-i-o.i 


27/32 


0.70 


3.33 


3.68 


2006 December 17 


1 


' oy -0.17 


_ „_un OQ 

u - u -o.o 


13.5/20 


0.86 


2.74 


3.11 


2007 January 01 


1 


✓M-.4-0 30 
i?u — 0.23 


^ ^ + 18 

u - u -o.o 


21.6/14 


0.09 


1.70 


2.27 


2007 February 13 


1 


|0 24 
■ y5 -0.19 


^ ^4-0 12 


17.5/14 


0.23 


1.90 


2.35 


2007 February 18 


1 


j-^ + 29 
■ ' y -0.25 


^ ^ + 13 

u - u -o.o 


16.5/12 


0.17 


1.76 


2.30 


2007 April 16 


1 


__-f 0.31 
' ' ' -0.21 


^ ^ + 15 

n n ' 
u - u -o.o 


10.9/15 


0.76 


2.02 


2.90 


2007 April 25 


1 


_ _|_0.32 
'^-0.27 


.-, .-,-1-0 20 

u - u -o.o 


13.7/16 


0.63 


1.92 


3.54 


2007 June 04 


1 


SQ"f~^ -^^ 
Oy -0.19 


1 ' ■ u - . 


10.7/16 


0.83 


2.31 


2.69 


2007 June 11 


2 


' 13 -0.27 


r-. .-,-1-0 2fi 




16.7/15 


0.34 


2.08 


1.97 


2007 July 24 


1 


on + 0. 28 
'° -0.18 


n n + 0.18 

O-0-o.o 


14.0/18 


0.73 


2.32 


3.22 


2007 July 30 


1 


^^-4-0 17 


^ „-|-0 07 

u - u -o.o 


16.8/24 


0.86 


3.28 


3.80 


2008 March 16 


1 


—^-1-0 09 
' ' y -0.08 


_ _-i_n 04 

u - u -o.o 


73.4/70 


0.37 


2.89 


4.05 


2008 April 13 


1 


■y°_o.o9 


,-,-t-O 03 

°-°Io.o 


58.5/73 


0.89 


3.26 


3.28 


2008 April 23 


1 


OO + 0.06 
b8 -0.09 


n n + 0.06 

°-°-o.o 


65.3/63 


0.40 


2.69 


3.24 


2008 June 01 


1 


' ' -0.10 


^ ^ + 04 
UiU -0.0 


55.5/45 


0.14 


1.78 


2.48 


20U8 July 05 


1 


nq + 0-11 


n n + 03 


b2.6/o2 


0.16 


2.32 


2.51 


2008 November 11 


1 


69+0 12 

u -0.10 


0+ 06 


42.4/46 


0.62 


1.75 


2.90 


2009 November 28 


1 


gg+0 .09 
uo -0.08 


0.03+° °, 6 


88.6/113 


0.96 


2.36 


4.28 


2010 February 09 


1 


7 o + 0-10 
'^-0.07 


o+ 07 


71.6/88 


0.89 


2.03 


3.20 


2010 April 17 


1 


gr+0.08 


o+ 05 


109.2/106 


0.40 


2.34 


4.14 


2010 June 25 


1 


71 +0-08 
' ' i -0.08 


o+ 03 


71.2/79 


0.72 


1.88 


2.96 


2010 November 11 


1 


gg+0.07 


o+ 04 


91.2/94 


0.56 


2.16 


3.57 


2011 January 21 


1 


DO -0.09 


+0 06 


54.6/54 


0.45 


1.18 


2.17 


2011 February 2 


1 


yo + 0.15 


0+ 07 


61.4/54 


0.23 


2.29 


3.35 



a All model fits include fixed, Galactic absorption of jVcai.H — 0.036 x 10 22 cm 2 (Dickey & Lockman 1990) and intrinsic (z — 0.658) absorption 7V H 

set as a free parameter. All errors arc at 90% confidence on one parameter. 

b P(x 2 1 v) is the probability of exceeding for u degrees of freedom if the model is correct. 
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TABLE 4 

Spectral Analysis of Image D of RX J1131-1231 



Observation 


r 


W H (z 3 =0.658) a 


X 2 /^ (cstat/i/) 




Soft Flux 


Hard Flux 


Date 




(10 22 cm -2 ) 




(10~ 13 erg 


£s cm~ 2 s _1 ) 


2004 April 12 


!-95Io 19 


o n 


(368/789) 




0.67 


0.73 


2006 March 10 




n n^ + 40 


(297/837) 




0.81 


1.13 


2006 March 15 


-■ ,-,1-1-1-0 49 


„ -1-0 32 
U - Z1 -0.21 


(290/837) 




0.78 


1.17 


2006 April 12 


- „ , |_ o 39 

1 - 65 -0 25 


f-, f-, + 26 

u - u -o.o 


(350/837) 




1.02 


1.46 


2006 November 10 


- .-^,-1-0 32 

!-52Io:i9 


^ ^ + 19 

O-OIo.o 


(403/837) 




1.11 


1.86 


2006 November 13 


O Q Q ~i~ » 5 3 
z ' dy -0.37 


^ ^ + 32 


(276/837) 




1.42 


0.59 


2006 December 17 


^ — 0.36 


,-,-1-0 26 

u - u -o.o 


(347.7/2008) 




0.83 


1.22 


2007 January 01 


- „ , i_ o 27 

A ' ao — 0.25 


f-, f-, + 06 

u - u -o.o 


(324.2/2003) 




0.99 


0.92 


2007 February 13 


2 - 10 -0.22 


r, ,-, + 0.18 

0-OIq.o 


(391.9/2009) 




1.45 


1.23 


2007 February 18 


_ +0.44 

2-01I .30 


^ ,-,+ 02 

0-0-oo 


(344.2/2009) 




1.05 


1.14 


2007 April 16 


„ „ __|_0 49 
Z - OO -0.42 


0-30-0.25 


(296.0/2009) 




0.79 


0.75 


2007 April 25 




„ „ + n ' > ] x 

u - u -o.o 


(329.3/2009) 




1.02 


1.27 


2007 June 04 


- ««4-n 40 
±- ' ^-0.34 


„ - ,,-1-0 27 

n 1 0^ 


(370.7/2003) 




0.80 


1.36 


2007 June 11 


„ - - -1-0 41 
Oil i u '^^ 
Z -- L1 -0.37 


_ _.__i_n 27 

0.243jl°- 2 I, 


(354.5/2009) 




0.96 


1.21 


2007 July 24 


„ __-l-0 40 
-0.27 


~ ~-l-0 22 
U ' U -0.0 


(313.4/2009) 




0.71 


1.05 


2007 July 30 


o n O + 0.47 
^• iM -0.26 


r, n + 0.20 

°-°-o.o 


(341.0/2009) 




1.00 


0.74 


2008 March 16 


_ _ „_|_0 1 6 

!-56Io is 


„ „-|-0 OR 
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a All model fits include fixed, Galactic absorption of -/Vcai.H — 0.036 x 10 22 cm 2 (Dickey 8z Lockman 1990) and intrinsic (z — 0.658) absorption 7V H 
set as a free parameter. All errors arc at 90% confidence on one parameter. 

b P(x 2 jv) is the probability of exceeding for v degrees of freedom if the model is correct. The C statistic docs not provide a goodness-of-fit measure. 
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Fig. 1. — Stacked image of all observations of RX J1131— 1231 listed in Table 1. The stacked image was filtered to include only photons 
with energies ranging between 0.4 and 3.0 keV and adaptively smoothed. There is no evidence for other sources in the field of view. The 
cross marks the center of the lens galaxy where we might see a fifth image. 
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Fig. 2. — The total (0.2-10 keV) light-curves of images A, B, C and D of RX J1131— 1231 shifted by the time delays estimated by Chantry 
et al. (2010). The total counts for images A and B have been corrected for pile-up effects and pile-up is unimportant for images C and D. 
The new X-ray data begin after epoch 6 (December 2006, JD — 2450000 = 4087). Data from epochs 1-6 were presented in Chartas et al. 
(2009) and Dai et al. (2010). 
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Fig. 3. — The total (0.2—10 keV, crosses with error bars) and SMARTS optical (i?-band, solid points) flux ratios A/B, C/B and D/B. 
The total band flux ratio C/B is relatively constant at 0.29 ± 0.03 supporting our hypothesis that the microlensing magnifications of images 
B and C during our observations are roughly constant. Bars in the top panel show the estimated optical (to) and X-ray (f-x) source crossing 
times for the size estimates from Dai et al. (2010) and the effective velocity estimate from Mosquera & Kochanck (2011a). The new optical 
and X-ray data begin after day 4053. The optical and X-ray data from days 2986 — 4053 were presented in Chartas et al. (2009) and Dai et 
al. (2010). 
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Fig. 4.— Soft (0.2-2 keV) and hard (2-10 keV) flux ratios D/B and D/C of RX J1131-1231. Significant energy dependent microlensing 
is detected in image D. 
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Fig. 5. — Simulated reflection spectrum from a cold accretion disk following Reynolds et al. (1998). The dashed line indicates the input 
incident power-law spectrum of the corona. We have over-plotted the rest-frame soft and hard bands used in our analysis. The soft band 
contains a larger fraction of direct emission than the hard band. The dominant fluorescence line is Fe Kct at 6.4 keV. 
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Fig. 6. — Stacked spectra of images C (panel a) and D (panels c) for the three periods shown in Figure 4. A close-up of the stacked spectra 
around the Fe Ka region for images C (panel b) and D (panel d). The curves are the best-fit absorbed power-law plus Gaussian line models. 
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Fig. 7. — Stacked spectra of image D shown in separate panels for the three periods shown in Figure 4. The curves are the best-fit absorbed 
power-law plus Gaussian line models. 
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Fig. 8. — The contours represent 68%, 90% and 99% x 2 confidence intervals of the flux normalizations, /p e Ka, of the detected Fe Ka line 
in image D (left panels) and in image C (right panels) in periods 1 (panels a and d), 2 (panels b and e) and 3 (panels c and f) from Figure 4. 
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Fig. 9. — The top panel shows the stacked spectrum of image C for all epochs fit by a power-law plus a Gaussian line model with Galactic 
absorption. The vertical dashed line indicates the energy of an Fe-Ka line with a rest-frame energy of 6.4 keV. The lower panel shows the 
residuals of the fit in units of la deviations. 
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Fig. 10. — Ratios of the (a) best-fit models (shown in Figure 6) and (b) observed flux densities D{(E)/C{E), where, D{(E) is the flux 
density for image D in period i and C(E) is the flux density for image C in all three periods. The three periods are shown in Figure 4. The 
model spectra are much sharper because they are not convolved to the instrumental resolution. 



